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Abstract 

The work studies wave activity in spatial systems, which exhibit non- 
local spatial interactions at the presence of a finite propagation speed. 
We find analytically propagation delay-induced wave instabilities for var- 
ious local excitatory and lateral inhibitory spatial interactions. The final 
numerical simulation confirms the analytical results. 

1 Introduction 

In the recent decades the space-time dynamics of extended complex systems has 
attracted much attention in physical, chemical and biological science [7|E1- I n 
this context the nonlocal interactions between system subunits have been stud- 
ied extensively 1161 . By virtue of this nonlocal spatial interaction, the 
propagation delay between two spatial locations may not be negligible if the the 
propagation speed in the system is small enough. Such propagation delays have 
been found to affect the activity spread solids an d plasmas jS], they af- 

fect the stability of neuronal systems or mav be responsible for subballistic 
transport in solids [22] • The present work addresses these phenomena and stud- 
ies the effect of a finite propagation speed in a generic spatio-temporal model. 
To describe the space-time evolution of nonlocal systems mathematically, partial 
differential equations have been studied vastly, e.g. Moreover integral- 

differential equations also have been studied in recent years to describe nonlocal 
interaction in spatially extended neuronal systems |21 151 1181 119U24) . Both latter 
equation types are strongly related for some specific integral-differential equa- 
tions [111 ^ . To illustrate this relation, assume a one-dimensional spatial field 
u(x) and the spatial kernel function K(x). Then it is 

/>oo 

K(x - y)u(y)dy = / K(k)u(k)e ikx dk 
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with the Fourier transforms K and u of K and u, respectively. Then choosing 
K{x) =exp(-\x\/VZ>)/2VD, it is K(k) = l/(l + Dk 2 ) = 1 -Dk 2 + D 2 k i + - ■ ■ 
and Eq. QJ reads 




K [x — y)u(y)dy s» i) + DAti(a;, t) 



+D 2 A 2 u(x, t) + ■ ■ ■ 

with A = d 2 /dx 2 . This means that fast-decaying integral kernels K(x) with 
V~D -C 1 represents a local diffusion process with diffusion constant D, while 
slowly-decaying integral kernels yield higher orders of spatial derivatives and 
thus represent long-ranged or nonlocal interactions. This procedure has been 
studied in some more detail in a previous work £Q and it turned out that integral- 
differential equations generalize some types of partial differential equations. In 
this context the present work extends the previous work by the study of a more 
general spatio-temporal model. 

In the research field of spatially extended systems, let the underlying model be 
partial differential equations or integral-differential equations, it is well-known 
that local excitation and lateral inhibition (local inhibition and lateral exci- 
tation) may yield stationary (non-stationary) spatial patterns in one or more 
spatial dimensions USES- However, a recent study showed that local inhibi- 
tion and lateral excitation may yield stationary instabilities as well for gamma- 
stributed spatial interactions 1 . This finding let us look more carefully on the 
nature of spatial interactions and the classification criteria for these instabil- 
ities. The present work investigates wether there are traveling waves in local 
excitatory and lateral inhibitory systems, which may contrast to the findings in 
previous studies. 

The work is structured as follows. The subsequent section introduces and mo- 
tivates the generic model studied and presents the steps of corresponding linear 
stability analysis about a stationar state. Then section applies the model 
to a field of harmonically damped oscillators, which are coupled by nonlocal 
excitation and inhibition. The analytically results are confirmed by numerical 
simulations. The subsequent summary closes the work. 



2 The generic model 

The model studied presumes identical elements which obey the (non-)linear 
differential equation 

f t V(x,t) = h[V(x,t)]. 

Here V(x,t) represents the scalar activity variable, T t = T t (d/dt) denotes the 
temporal operator and the functional h[V(x, t)] defines the local dynamics at 
spatial location x. Let the elements {x} be located in a spatial domain and 
coupled to each other. This configuration reflects a topological network with 
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identical elements in which no spatial metric is considered, or a topological net- 
work embedded in a physical space [211 OJ. In the the present model the single 
elements are connected according to a probability density function similar to 
topological networks but their connectivity function depend on the Euclidean 
distance between the elements similar to spatial networks. Further the spatial 
coupling is presumed homogeneous in space, that is the coupling is dependent on 
the spatial distance between the elements. These properties describe a scale-free 
lattice Moreover, the presented model presumes infinitely dense ele- 

ments, that is we consider a continuous field. This contrasts to networks whose 
major element is the discreteness. However in biological systems, e.g. in the 
brain, the single cells are dense and the continuum limit holds . Furthermore 
the speed by which the activity is transmitted between single elements shall be 
finite. This finite propagation speed also contrasts to most spatial networks, 
which frequently neglect propagation delay effects. 

Considering all previous aspects, the clement at spatial location x obeys the 
integral-differential equation 



Here I(x, t) represents the external input and the kernel functions K, L rep- 
resent the probability density functions of spatial connections in two networks, 
and thus K{x) and L(x) are normalized to unity in the spatial domain J7. 
In the following we assume fi = TZ. The functionals f[V] and g[V] repre- 
sent the (non-)linear coupling functions. Further the model considers prop- 
agation delays between spatial locations x and y at distance \x — y\ due to 
the finite transmission speed vk and vl corresponding to the two networks. 
This model allows the study of a wide range of combined spatial interaction. 
For instance spatial systems involving both mean-field and local interactions 
may be modeled by K(x) = const and L(x) = S(x), respectively. The sub- 
sequent paragraphs focus on the combination of excitatory and inhibitory in- 
teractions with kernel functions K(x) and L(x), respectively. In addition we 
specify f[V] — cikS{V], g[V] — o,lS{V] with the nonlinear functional S[V] and 
the total amount of excitation and inhibition <ik and a l , respectively. Further 
we presume cik > cll- For h(V) = 0, Eq. represents a well-studied nonlocal 
model of neural populations fM l3l ITO! IT1 15] . 

First let us determine the stationary state which is constant in space and 
time. Applying a constant external input 1$, we find the implicit equation 
= h(Vo) + f(Vo) + g(Vo) + Iq. Here we used the fact that the kernel functions 
are normalized to unity. Then small deviations u{x, t) about the stationary 
state obey a linear evolution equation and relax in time to the stationary state 
V with u(x,t) ~ e At , A e C if Re(X) < 0. Then the subsequent expansion of 



T t V{x,t) 



h[V(x,t)] + I(x,t) 





+L(x-y)g V(y,t 
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u(x,t) in the continuous Fourier basis {e lkz } yields the implicit condition 

T(A) = h Q + S' I dy 
Jn 

x {a K K{y)e- x ^' VK - a L L(y)e- x ^/ VL ^j e lky 

(1) 

with ho = dh/dV, S' — dS/dV computed at V = Vq. Since the external 
input drives the system, its stability is subjected to Iq and subsequently 5" 
represents the reasonable control parameter. Equation is difficult to solve 
exactly for A. In a previous study a more specific model has been studied in 
detail for stationary and non-stationary instabilities pQ . We follow its approach 
and consider large but finite propagation speeds. Then it is exp(— X\z\/v) « 
I - (X/v)\z\ for \X\/vjc, \X L \/v L -C I and reads 

T(X) ^h + S'Ko(k) - XS'K^k) . (2) 

with 

K (k) = a K K^{k)-a L m{k) 



with fWffc) = F(z)\z\ n exp(—ikz)dz. The term Ko{k) and K\(k) repre- 
sents the Fourier transform and the first kernel Fourier moment of cikK(x) + 
(IlL(x), respectively. Essentially let us specify the temporal evolution of a single 
element to a damped oscillator with 

t d 2 d 

Then the stability threshold of stationary and non-stationary instabilities read 
1 K (k c ) 



S' c l -ho 

i_ _ k^k) 

S'c 7 



(stationary instability) (3) 
(non-stationary instability) 

(4) 



Hence the global maximum of K (k c ) and — K\(k c ) define the critical wavenum- 
ber k c and the critical threshold S' c of stationary and non-stationary instabilities, 
respectively. Figure [T] illustrates this result for non-stationary instabilities. In 
the following, we call an emerging non-stationary phenomena a wave instability, 
if k c ^ 0. and the instability with k c = is called a global oscillation. Moreover 
it is interesting to note that the threshold of the non-stationary instability does 
not depend on the local (non-)linear dynamics of the elements h(V). 
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Figure 1: Sketch to illustrate the mechanism of non-stationary instabilities. The 
plot shows the right hand side of Eq. (@J (solid line) and its left hand side (dot- 
dashed line). If the stability threshold is not reached, i.e. S' = s< < S' c , the 
system is stable. In contrast, if S' — s c = S' c (S' — s> > S' c ) the system is 
marginal stable (unstable). 



3 Application to specific model 

First let us discuss briefly various types of nonlocal interactions. We choose 
the excitatory and inhibitory spatial kernel to the gamma function and the 
decreasing exponential function, resp., with 

K(z) = - * , \z\ p - 1 e-^ q , L(z) = ^-e~W r (5) 
v ; 2qPT(p) 1 1 ' w 2r w 

Since the Eq. Q does not depend on h[V], we set h[V] = 0. In other words, the 
single elements in the field are assumed damped harmonic oscillators, which are 
coupled by an excitatory and an inhibitory network. The parameter r defines 
the spatial range of the inhibitory interaction, while both p and q represent the 
spatial constants of the excitatory nonlocal interaction. The scase p = 1 is well- 
studied, see e.g. |1(J) . and the combination of q and r yield the four important 
cases, namely the global excitation, the global inhibition, the local excitation- 
lateral inhibition and the local inhibition-lateral excitation. Further the choice 
p > 1 yields local inhibitory interactions or local inhibitory-lateral excitatory 
interactions pQ. The present work focus on the case p < 1 which yields the 
local excitation or the local excitation-lateral inhibition. Please see Fig. [5] for 
the corresponding illustration of the cases p < 1 and p > 1. In order to gain 
the critical wavenumbcr of the wave instability, the kernel Fourier moment in 
Eq. @ is computed to 

_ pqcos((p+ l)arctan(fcg)) 
{) ~ (l + gajfe2)(H-D/2 W 

_ rcos(2arctan(fcr)) 
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Figure 2: Illustration for some spatial interaction types for different parameters 
p, q and r. Panel (a) shows the kernel function for p < 1, q > r that yields global 
excitations while panel (b) reveals local excitation-lateral inhibition for p < 
1, q < r. In panel (c) the kernel function has been illustrated for p > 1, q < r, 
which represents global inhibition, while panel (d) shows local inhibition-lateral 
excitation for p > 1, q > r. It is M(x) = K(x) — L(x) 

Since K±(k) — > for |fc| — > oo, there is a maximum of — at some |fc| = |fc c | 

with -Kx(k c ) > if -i^i(O) > and -(fk^/dk 2 > 0. These sufficient 
conditions read 

?>pK,( 3 <^p(p+1)(p + 2) (8) 
ba r 

with the new parameters £ = r/g, w r = Vl/vk and a r = a^/a^ < 1. Figure|31 
illustrates the different parameter regimes of the conditions. At first it turns 
out that there is large parameter regime which allows for non-stationary insta- 
bilities. In order to gain further insight to the nature of spatial interactions for 
p < 1, Fig. 0] shows the parameter regimes and the different number of roots 
of the kernel function K(z). We find a parameter regime of single roots, which 
reflects local excitation-lateral inhibition interaction, and a parameter regime of 
no roots, i.e. global excitation. Moreover, there is a parameter regime of mul- 
tiple roots which, for instance, reflects local excitation - mid-range inhibition - 
lateral excitation and allows for wave instabilities. 

To illustrate the relation of the spatial interaction type and the resulting 
non-stationary instability, Figure [S] plots the kernel functions K{z) and the 
corresponding negative first kernel Fourier moments —K\{k) for four different 
cases of spatial interactions. These specific cases are denoted by filled dots in 
Fig. 01 We find local excitation-lateral inhibition for £ = 0.8 (top row, left panel 
in Fig. |2J|. Since the global maximum of — K\(k) gives the critical wavenum- 
ber k c , this case yields global oscillations with k c = (top row, right panel 
in Fig. |3J|. For £ = 0.75, now the kernel function exhibits two roots which 
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Figure 3: The occurrence of oscillatory instabilities subject to the parameters 
£ = r/q and p. The dashed and thin solid line represent the function £ = 
(v r p(j)+ 2)/6ct,-) 1 / 3 and £= p/a r , respectively, taken from Eq. JHJ). Here 

it is v r = 1, a r = 0.8 . 




Figure 4: Number of zero-croosings, i.e. roots, of the kernel function K{z) and 
the corresponding nature of spatial interactions subject to the parameters £ and 
p. The dashed and thin solid line are taken from Fig. [3J The red lines delimits 
the novel parameter regime of multiple roots. In addition to the parameters 
from Fig. it is clk = 10, — 8, vk = vl = 1 and q = 1. 
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Figure 5: The kernel functions axK(z) + a^L(z) and the negative first kernel 
Fourier moment K\{k) for the four different cases taken from Fig.|H The insets 
in the left coloumn represent foci of the corresponding plot to illustrate the 
additional lateral excitation. The parameters are taken from Fig. 0] and the 
different values of £ correspond to the different dots in Fig. 0] 

still leads to global oscillations. In case of £ = 0.70 the sufficient condition for 
wave instabilities are fulfilled, as —Ki(k) reveals a critical wavenumber k c 0, 
see Fig. 0] Moreover, the spatial interaction function reveals local excitation, 
mid-range inhibition and lateral excitation, which has not been found yet in 
previous studies. At last for £ = 0.55 there is also a wave instability showing 
local excitation, mid-range inhibition and lateral excitation. In addition to this 
study, Fig. compares the spatial interaction function, the Fourier transform 
K and — K\{k) for two values of £. This comparison is necessary as the sta- 
tionary instability may occur first before the wave instability while increasing 
the control parameter from small values. For both parameter values of £ the 
global maximum of — K\{k) exceeds the global maximum of Kq and thus a wave 
instability and not a stationary instability may occur. Thus Fig. supports the 
previous findings. 

To visualize and confirm the analytical findings, finally we integrate the 
integral-differential equation Q numerically. The Euler integration scheme 
is applied for the time evolution integration and the Monte-Carlo method VE- 
GAS ^3 E] computed the spatial integration. This specific Monte-Carlo inte- 
gration method combines stratified and importance sampling, which is impor- 
tant to gain good estimates of the divergent excitatory kernel function. Moreover 
boundary conditions are set periodic with the period L while the initial condi- 
tions have to be defined in the time interval to G [— L/v;0]. We have chosen 
V(z,to) = cos(k c x) + 0.5 * cos(5fc c x). Figure shows the wave instability for 
p = 0.5, £ = 0.7 discussed above, cf. Fig.Q] and confirms the analytical finding. 
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Figure 6: The comparison of the kernel functions K(z) + cllL(z) (left 
coloumn), the Fourier transform Ko(k) (dashed line in right coloumn) and 
— K\(k) (solid line in right coloumn) for two values of £. Additional param- 
eters are taken from Fig. 




Figure 7: Space-time plot of the integrated activity for p = 0.5, r = 0.7 shows 
the wave instability for local excitation - lateral inhibition interactions. Ad- 
ditional parameters are cik — 10, ah = 8, vk = vl = 0.6, 7 = 0.2, Io — 
2.8, dx = 0.089, dt = 0.08 and h[V] = 0, S[V] = 1/(1 + exp(-1.84(V - 3.0))). 
These parameters yields k c — 0.4. 
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4 Summary 



The present work describes the nonlocal spatial interaction of single elements 
in a continuous field. The application to damped oscillators, which are coupled 
by both nonlocal excitation and inhibition, reveals a new mechanism of wave 
instabilities. In contrast to previous findings in pattern forming systems, we 
find wave instabilities for local excitation and lateral inhibition. To obtain more 
information on the corresponding spatial interaction, we investigated addition- 
ally the spatial interaction function. It turns out that various spatial interaction 
types exhibit wave instabilities, as local excitation-lateral inhibition and local 
excitation-midrange inhibition-lateral excitation. Hence the sign of the spatial 
interaction function does not carry the full information of the expected instabil- 
ity patterns in spatial systems with finite propagation delay. We have observed 
that the comparison of the global maximum of the Fourier component and global 
maximum of the negative first kernel Fourier moment allows for a reliable esti- 
mation of the instability type. 

Future work shall study the classifiers in vectorial models, as the nonlocal version 
of the Bruesselator, and the classification scheme in two- and three-dimensional 
systems. 



Appendix: The kernel Fourier moments 

The n-th kernel Fourier moment of the specific kernels JSJ can be computed 
from 



/>oo 

J(p) = / x p - 1 e- z/q cos(kx)dx 
Jo 



1 \ ~ p/2 

= T(p) ( — + k 2 \ cos (parctan(fccjr)) 



q 2 

Then we find for the excitatory kernel 

qPT(p) 

T(p + n) q n cos ((p + n) arctan(fc<7)) 

T(p) ( 1 + g 2 fc 2)(p+n)/2 • 

In addition, we find 

-^K {n \k) = k i - n+2 \k) . 

Since the inhibitory kernel is a specific case of the excitatory kernel with p = 1 
and q = r, we obtain Eqs. ©, J7J) from the previous expressions. 
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